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ABSTRACT 

We observed the young energetic pulsar J1617-5055 with the Chandra ACIS detector for 60 ks. In addition 
to the pulsar, the X-ray images show a faint pulsar- wind nebula (PWN) seen up to ^ 1' from the pulsar. De- 
convolution and reconstruction of the image reveal a brighter compact PWN component of ^ 1" size, possibly 
with a jet-torus morphology. The PWN spectrum fits an absorbed power-law (PL) model with the photon index 
r « 1.5. The total PWN luminosity in the 0.5-8 keV band, Lp„n = 3.2 x 10^^ ergs s"' for d = 6.5 kpc, is a frac- 
tion of 2 X 10 of the pulsar's spin-down power E and a fraction of 0.2 of the pulsar's X-ray luminosity, which 
is a factor of 20 lower than one would expect from an average empirical relation found for a sample of PWNe 
observed with Chandra. The pulsar's spectrum can be described by an absorbed PL with nn » 3.5 x 10^^ cm~-^ 
and r « 1 . 1 , harder than any other pulsar spectrum reliably measured in the soft X-ray band. This non-thermal 
emission is » 50% pulsed showing one peak per period. We have also investigated a possible connection be- 
tween J1617 and the extended TeV source HESS J1616-508 whose center is located 10' west of the pulsar. We 
find no preferential extension of the X-ray PWN toward the TeV source. Therefore, the Chandra data do not 
provide conclusive evidence for PSR J1617-5055 and HESS J1616-508 association. We have also analyzed 
archival X-ray, radio, and IR data on the HESS J1616-508 region and found traces of diffuse emission (re- 
sembling a shell in the radio) coinciding with the central part of HESS J1616-508. We speculate that the TeV 
source may be multiple, with most of the emission coming from an unknown SNR or a star-forming region, 
while some fraction of the TeV emission slill may be attributed to the J1617 PWN. 

Subject headings: ISM: individual (HESS J1616-508) — pulsars: individual (PSR J1617-5055) — stars: 
neutron — X-rays: ISM 



1. INTRODUCTION 

High-resolution X-ray observations of young (r < 10 kyrs) 
rotation-powered pulsars usually show a pointlike pulsar em- 
bedded in a pulsar wind nebula (PWN). In most cases, the pul- 
sar emission is dominated by the nonthermal emission gener- 
ated in the pulsar magnetosphere. The X-ray spectrum of this 
emission usually fits power-law model with typical photon in- 
dices r ~ 1-2. As a rule, the nonthermal emission component 
is strongly pulsed, showing one or more peaks per period. 
Studying the non-thermal emission constrains the emission 
and particle acceleration mechanisms in the pulsar magneto- 
sphere. Young pulsars are also interesting because they emit 
strong magnetized relativistic winds which interact with the 
ambient medium and form PWNe whose emission is inter- 
preted as synchrotron radiation from the shocked pulsar wind 
(see Kaspi et al. 2006; Gaensler & Slane 2006; Kargaltsev & 
Pavlov 2008, hereafter KP08, for reviews). Studying PWNe 
helps to understand the structure and dynamics of relativistic 
pulsar winds, elucidate the mechanisms of PWN formation, 
evolution and interaction with the ambient medium, and es- 
tablish the properties of the relativistic plasma in PWNe. 

The young (r = P/2P = 8.13 kyr; £ = 1.6 x lO^'^ ergs s'^) 
69 ms pulsar J1617-5055 (hereafter J1617) was first identified 
through its X-ray pulsations with ASCA (Torii et al. 1998)-^ 
with the radio pulsations found shortly afterwards (Kaspi et 

' Dept. of Astronomy, University of Florida, Bryant Space Science Cen- 
ter, Gainesville, FL 3261 1; oyklOO@astro.ufl.edu 

^ The Pennsylvania State University, 525 Davey Lab, University Park, PA 
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' The 69 ms pulsations from this region were first detected with Ginga 
(Aoki et al. 1992), but it was impossible to determine the source of these 
pulsations because of the very poor angular resolution of Ginga. 



al. 1998). The ASCA spectrum of the pulsar was described 
by the absorbed power-law (PL) model with the photon index 
r = 1 .4 ± 0.2 and observed flux Fx = (3 .6 ± 0.2) x lO^^^ ergs 
s-i cm-2 in the 3-10 keV band (Torii et al. 2000). 

Since the ASCA detection, J1617 has been observed with 
other X-ray satellites. It was imaged 6' off-axis in Chan- 
dra ACIS observations of the RCW 103 supernova remnant 
(SNR) in 1999 September (Garmire et al. 1999) and Febru- 
ary 2000. It was also observed by XMM-Newton in 2001 
September (two pointings, 30 ks aimed at the pulsar, and 20 
ks centered at RCW 103). The J 1617 properties inferred from 
these XMM-Newton observations have been briefly reported 
by Becker & Aschenbach (2002; hereafter BA02). They 
found that the spectrum in the 0.5-10 keV band fits the ab- 
sorbed PL model wifli F = 1.1-1.4, Hh = (2.8-3.6) x lO^^ 
cm"^, and unabsorbed flux = (4,9_5.4) x 10^'^ ergs"^ 

cm"^ (the numbers correspond to the 90% confidence range). 
BA02 have also detected pulsations with a single, asymmetric 
pulse and a pulsed fraction of ^ 50% in the 2.5-15 keV band. 
These observations confirmed the non-thermal nature of the 
pulsar emission; however, they lacked the resolution needed 
to separate the pulsar emission from the emission of a possi- 
ble compact PWN around this young, remote pulsar J 1617 
was also observed serendipitously by XMM-Newton in 2005 
August (w 90 ks). 

J 16 17 has been also detected at higher energies with RXTE 
PCA/HEXTE (2-60 keV band; Torii et al. 2000; Kuiper 
2007), and Integral IBIS/ISGRl (20-300 keV band; Landi 
et al. 2007). Landi et al. (2007) have produced a joined 
spectral fit to the XMM-Newton, BeppoSAX, and Integral 
data and found fitting parameters close to those obtained 
by BA02 from the XMM-Newton data alone [F = 1.4± 0.1, 
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hh = (3.9 ±0.3) X 10^2 and F™*^ = 4.2 x IQ-'^ ergs 

cm"^ in 2-10 keV]. Fitting the RXTE spectrum in the 
2-30 keV band with a power-law model, Kuiper (2007) found 
r= 1.30 ±0.01 for fixed nH = 3.2 X 10^^ cm'^. 

Observations of J1617 and its surroundings have become 
particularly interesting after the discovery of the extended 
(- 20' in diameter) TeV source HESS J1616-508 (hereafter 
HESS J1616) by Aharonian et al. (2006). Although its center 
of gravity is offset by w 10' - 1 1' from the pulsar, the extent of 
the TeV emission still encompasses the pulsar location. The 
HESS J 16 16 spectrum above 200 GeV fits a power-law model 
with r = 2.35 ±0.06 and 1-10 TeV flux Frev w 1.7 x lO"" 
ergs s"' cm"^. The central region of HESS J1616 was ob- 
served with Suzaku XIS, but no X-ray counterpart was found 
to a limiting flux of 3 . 1 x lO"'"* ergs cm"^ s"' , so that the TeV 
source was suggested to be "the best example in the Galaxy 
of a dark particle accelerator" (Matsumoto et al. 2007). Landi 
et al. (2007) observed the field of HESSJ1616 with Inte- 
gral IBIS/ISGRI and Swift XRT, and also analyzed archival 
XMM-Newton and BeppoSAX data. These authors con- 
cluded that J1617 is the only plausible counterpart to the TeV 
source. 

For many young pulsars, radio or X-ray SNR associations 
have been found, but no SNR association has been established 
for J1617 so far. Although J1617 is just ^ 7' north of the cen- 
ter of the RCW 103 SNR, outside of the 5' SNR radius, the 
pulsar is almost certainly unrelated to the SNR. RCW 103 
hosts a central compact object (CCO; possibly a NS binary or 
a transient magnetar; e.g., Pavlov et al. 2002, 2004; de Luca et 
al. 2007), which is a more plausible candidate for the compact 
remnant of the SN explosion. In addition, the RCW 103 spec- 
trum shows a significantly smaller interstellar absorbtion than 
J1617 (Gotthelf et al. 1997). The relative positions and sizes 
of HESS J1616 and RCW 103 suggest that the HESS source 
is also not associated with the SNR. 

To search for a compact X-ray PWN around J1617 and sep- 
arate the pulsar and PWN emission, we have carried out a 
deep Chandra ACIS observation of J 1617 in 2006 June, with 
the pulsar placed near the optical axis. Here we report the 
results of this observation and our analysis of archival mul- 
tiwavelength data on J1617 and HESS J1616. In §2 we de- 
scribe the observations, the images of the J1617 PWN and the 
HESS J1616 field, and the spectral analysis of the pulsar and 
the PWN. We discuss the pulsar and PWN properties and the 
nature of HESS J1616 in §3, and summarize our findings in 
§4. 

2. OBSERVATIONS AND DATA ANALYSIS 

We observed J 16 17 with the Chandra ACIS detector on 
2006 June 6 (ObslD 6684) for 60 ks in Timed Exposure 
(TE) mode with Very Faint (VF) telemtry format. The pul- 
sar was imaged on the AClS-13 chip near the optical axis. 
The other ACIS chips were turned off. To reduce the pile-up 
in the pulsar image, we used 1/4 subarray 8' x 2' field of 
view [FOV]; frame time 0.84104 s, including 0.8 s exposure 
time and 0.04104 s "dead time"). No substantial background 
changes were detected during the entire observation. The sci- 
entific exposure time was 57.24 ks. 

We will also use in this work two archival Chandra ACIS 
imaging observations of 1999 September 26 (ObslD 123) and 
2000 February 8 (970). As the primary target of those obser- 
vations was the RCW 103 SNR, the pulsar was imaged w 6' 
off-axis, which resulted in a broadened point spread function 
(PSF). The frame time was 3.24104 s in both observations 
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Fig. 1. — Chandra images of the field around J1617 with the observation 
dates listed in the images (see Table 1 for details). The bottom panel shows 
the combined image from 1999 and 2000 observations. All the images ai'e 
binned to a pixel size of 4" and then smoothed with the r = 12" gaussian 
kernel. The contrast is chosen to emphasize the faint large-scale emission. 
The extended shell, partly seen in the middle and bottom images, belongs to 
RCW 103. 

(other details are given in Table 1.) The second observation 
suffered from significant background flares (up to 6 times the 
quiescent level). Filtering out the flares reduced the exposure 
by w 5 ks. 

In addition we analyzed the ACIS observation of of 2002 
March 3 (ObslD 2759). The data were taken in Continuous 
Clocking (CC) mode, which allows one to achieve time reso- 
lution of 2.85 ms at the expense of spatial information in one 
dimension. 

The central part of HESS 1616 was serendipitously imaged 
on the ACIS S2 and S3 chips during the Chandra observations 
of the Kes 32 SNR that occurred on 2001 October 10 (ObslD 
1960). The data were taken in Faint Mode, and the scientific 
exposure time was 29 ks (Vink 2004). 

Reduction of the Chandra data was done with the Chandra 
Interactive Analysis Observations (ClAO) software (ver 3.4, 
4.0; CALDB ver 3.3.0.1, 3.4.0) and FTOOLS (ver. 6.3). In 
our image analysis we have made use of MARX'* and Chandra 
Ray Tracer (ChaRT) software''. We also used XSPEC (ver. 
1 1.3.2) for the spectral analysis. 

J 16 17 was first observed with the XMM-Newton EPIC 
MOSl and MOS2 detectors on 2001 September 9 (ObslD 

MARX (Model of AXAF Response to X-rays) is a suite of programs de- 
signed to enable the user to simulate the on-orbit performance of the Chandra 
satellite. See http://space.mit.edu/ASC/MARX/ 

^ The sofware is available at http://cxc.harvard.edu/chai1/. 
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TABLE 1 
Chandra Observations of J1617. 



ObsID 123 970 2759 6684 

Date 1999 Sep 26 2000 Feb 8 2002 March 03 2006 June 6 

Mode, chip region TE, Full Frame TE, Full Frame CC, Full Frame TE, 1/4-subarray 

Telemetry format VF F F VF 

Chips activated 10-13, S2,S3 I2,I3,S2-S4 I2,13,S1-S4 13 

J1617 imaged on 11 13 S2 13 

OiT-axis angle 6.1' 6.1' 6.3' 40" 

Sci. Exp.,ks 13.92 13.75 48.82 57.24 



Note. — Except for ObsID 970, the scientific exposures coincide with the live times given in the EVT2 file headers (for the chip of interest). In ObsID 970, 
removing the flares reduced this time by 5 ks. 




TABLE 2 

Count statistics for the spectral extraction regions 



Region 


A 


Moi 






S/N 


S 


inner PWN 
outer PWN 


3.14 
2610 


1162 ±32 
959±31 


732 ±75 
445 ±21 


430 ±82 
514±37 


5.2 
13.9 


2.39 ±0.46 
(3.44 ± 0.25) X 10-3 



Note. — Source (Nsm), total (Ntod and background (Wbg) counts are ex- 
tracted from the regions with an area A (in arcsec^) shown in Figures 2 and 
4, in the 0.5-8 keV band. For inner PWN the uncertainty of the background, 
~ 10%, is mostly systematic (due to the uncertainties involved in the PSF 
simulation). The mean surface brightness cS is in units of counts ks"' arcsec"^ 
in the 0.5-8 keV band. 

0113050701)^. J1617 was also serendipitously observed with 
the same detectors during the observation of RCW 103 on 
2005 August 23 and 24 (ObsID 0302390101). In both ob- 
servations the two MOS cameras were operated in Full Win- 
dow mode (2.6 s time resolution) with the Medium optical 
filter; the dead-time corrected exposure times were w 28 and 
K, 86 ks. During the 2001 observation, the EPIC PN camera 
was in Timing mode, while it was in Small Window mode, 
with J 1617 outside of the FOV, during the 2005 observation. 
For the XMM-Newton data analysis, we used the Processing 
Pipeline Subsystem (PPS). Both the 2001 and 2005 observa- 
tions were strongly contaminated by flares during which the 
background count rate exceeded the quiescent level by a fac- 
tor of up to 16. Filtering the flares out left 22.2 ks and 70.3 ks 
useful scientific exposure for the first and second observations 
respectively. 

2. 1 . PWN and pulsar images 

Figure 1 shows binned and smoothed large-scale images of 
J1617 and its surroundings in the 2-8 keV band. The upper 
panel shows the 1/4 suban^ay AClS-13 image from our obser- 
vation of 2006 June. In this image we see a relatively bright 
pointlike source at the radio pulsar position embedded in faint 
diffuse emission more extended toward south and southeast of 
the pulsar The asymmetry is better seen in the ACIS images 
from earlier observations (Fig. 1, middle and bottom), which 
have larger FOVs, with the pulsar w 3' north of the RCW 103 
shell seen near the bottom in each of the two images. 

Figure 2 shows a close-up view of J1617 and its immediate 
vicinity from our observation of 2006 June. With the pulsar 
imaged on-axis, this observation provides a much better reso- 
lution than the off-axis ACIS observations of 1999 and 2000. 
The pointlike source, surrounded by diffuse emission, is cen- 

' J1617 was also serendipitously imaged off-axis in the observation 
0113050601 on the same date. We do not analyze this observation because 
a much longer XMM-Newton observation of 2005 August is available, with 
the about the same pointing and roll angle. 



Fig. 2. — ACIS-I3 (1/4 subarray) image of J1617 and its vicinity in the 
2-8 keV band, binned by a factor of 2 (pixel size 0.98"). Top: Unsmoothed 
image with the regions used for extraction of the extended PWN spectrum 
(within the polygon, excluding the circle; see §2.3.2 and Tables 2 and 3) and 
the background (box). Bottom: The same image adaptively smoothed with 
csmooth. 



tered atR.A.= 16*'17"29.35% decl.= -50°55'12.8". The posi- 
tion uncertainty is dominated by the uncertainty of the Chan- 
dra absolute astrometry, 0.65" at the 90% confidence level 
for sources within 2' from the optical axis on the ACIS -13 
chip^. As this position differs by only 0.62" from the radio 
position of J1617 reported by Kaspi et al. (1999), the point- 
like source must be the X-ray counterpart of the radio pul- 
sar. The amorphous diffuse emission around J1617 is dis- 
tinguishable from the background up to w 1' from the pul- 
sar. The background-subtracted surface brightness of this 
extended PWN is ~ 0.002-0.01 counts arcsec"^ ks"'. 



being 



' See §5.4 and Fig. 5.1 of the CAandra Proposers' Observatory Guide v. 10 
at http://asc.harvard.edu/proposer/POG. 
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Fig. 3. — Radial profiles of the observed emission around the pulsar po- 
sition and the simulated PSF in the 2-8 keV band. The histogram in the top 
panel (bin size 1") demonstrates the large-scale PWN emission. The his- 
togram in the bottom panel (bin size 0.25") shows the compact inner PWN 
contribution, including the excess at 0.5"-1.5" from the pulsar (see also 
§2.1 and Fig. 4). 

somewhat brighter south and southwest of the pulsar. 

The radial profile of the detected emission, centered on the 
pulsar position, is shown in Figure 3, together with the sim- 
ulated PSF^ and backgroimd level. We see that the extended, 
faint PWN prevails over the pulsar emission at distances > 5" 
from the pulsar. 

Given the large distance to the pulsar and the faintness of 
the extended PWN, one could expect a more compact PWN 
component, which would not be immediately seen in pipeline- 
processed data. To search for such a component, we produced 
images of the immediate pulsar vicinity at subpixel resolu- 
tion (see Fig. 4). We first removed the pipeline pixel random- 
ization and appUed the sub-pixel resolution tool to split-pixel 
events in the image (Tsunemi et al. 2001; Mori et al. 2001). 
We then simulated the PSF image in the 2-8 keV energy band 
at the same position on the detector as in the real observation. 
We binned the observed image and the simulated PSF by a 
factor of 0.2 compared to the original ACIS pixel size. Fi- 
nally, we used the CIAO ver. 4.0 implementation (provided as 
anarestore script) of the Lucy-Richardson deconvolution 
algorithm (Lucy 1974). We ran arestore with several set- 
tings for different numbers of iterations, between 30 and 200, 
and noticed that the iterations converge (i.e., no new structure 

^ To simulate the PSF with ChaRT and MARX, we followed 
the steps outlined in http://cxc.harvard.edu/chart/threads/index.html. We 
have also tried several values of the Dither Blur parameter (see 
http://cxc.harvard.edu/chart/caveats.html and Misanovic et al. 2007 for dis- 
cussion) and adopted the value of 0.1" that provides the best match between 
the simulated PSF and the real data within r < 3" from the pulsar. 



appear) once the number of iterations approaches 100. 

Examples of such deconvolved images are shown in Figure 
4 for 30, 100, and 200 iterations (panels b, c, and d, respec- 
tively). When the number of iterations approaches 30, the de- 
convolved image starts to reveal an extended structure (~ 2" 
in linear extent) elongated in the south-north direction, with 
the point source being in the middle of this structure. With in- 
creasing the number of iterations, the elongated structure be- 
comes one-sided (i.e., it appears only to the south of the point 
source), and an arclike feature, roughly perpendicular to the 
southern elongation, appears at about 1" north of the pulsar. 
Increasing the number of iterations beyond w 100 does not 
change the structure of the compact extended emission. The 
exposed features of the deconvloved images could be inter- 
preted as a jet (south of the pulsar) and an arc (perhaps part 
of a ring or a torus). Note that similar features are often seen 
in Chandra images of PWNe around young pulsars (see e.g., 
KP08). The radial profile shown in Figure 3 (bottom) demon- 
strates the significance of the count excess at r w 0.5"-1.5" 
from the pulsar. 

For comparison, we also show a deconvolved ACIS image 
without removing the pipeline pixel randomization (Fig. 4e). 
Although the deconvolution of this image does not recover 
as much structure as the images with the pixel randomization 
removed, it does show some of the features that are seen in 
the deconvolved non-randomized images. 

To test reliability of the deconvolution procedure, we ap- 
phed it to the ACIS image of the CCO in the Cas A SNR (Ob- 
slD 6690), a radio-silent NS showing no radio pulsar activity, 
for which no PWN is expected (e.g., Pavlov et al. 2004). The 
ACIS image of the CCO has « 7200 counts within the 0.98" 
radius around the source, comparable to J1617 in our ACIS- 
1 image. Unlike J 1617, the deconvolved image of the Cas 
A CCO preserves the pointlike appearance with no extended 
structure (see Fig. 4f). Although, as any deconvolution al- 
gorithm, arestore could introduce artifacts in the images, 
the test example of the Cas A CCO provides credibility to the 
features seen in the deconvolved images of J 1617. 

2.2. The field of HESS J 1616-508. 

Since J1617 is considered as a plausible counterpart for 
the extended TeV source HESS J1616 (Landi et al. 2007), we 
have searched for additional data in the Chandra and XMM- 
Newton archives that would provide at least partial coverage 
of the HESS J 16 16 extent and allow us to look for other po- 
tential X-ray counterparts of the TeV source and examine its 
connection to J1617. In the Chandra archive we found an 
ACIS observation of the Kes 32 SNR (Vink 2004). During 
that observation, the S2 and S3 chips imaged the central part 
of HESS J1616 (see Fig. 5a). The brightest source (marked X 
in Fig. 5a) is located on the S2 chip around R.A.« 16'*16'"10^ 
decl.ss -50°54'30". Because of the small number of counts 
(140 ±20 counts within an r = 40" circle in 0.5-8 keV; 
S/N « 7) and the broad PSF (FWHM w 15" at the off-axis 
angle 9 ~ 16'), it is difficult to determine whether the source is 
extended or multiple; however, with the linear extent of ~ 1', 
it does not appear to be point-like. Assuming an absorbed 
PL spectrum with F = 1.5 and nu,22 = ne/lO^^cm"^ = 3.45 
(as found from the fit to the pulsar spectrum; see §2.3.1), the 
unabsorbed 0.5-8 keV flux of source X can be estimated as 
Fx ~ (1.7-2.2) X 10-13 ergs cm'^ s'K 

We see some hints of extended emission southwest of 
source X; however, the low count rate, ~ 1.4 counts ks"' 
arcmin-^ in 0.5-8 keV (corresponding to the unabsorbed sur- 



5 




J1617 



Fig. 4. — a: Original ACIS-I3 image of pulsar and its immediate vicinty (pixel size 0.1"), with pipeline pixel randomization removed (circle radii are 0.5", 
0.74" and 1.25"). b, c, and d: Deconvolved images obtained with arestore from the image in panel a after 30, 100, and 200 iterations, respectively, e: Same 
as c but without removing the pipehne pixel randomization in the original image. /■ Original (top left; pixel size 0.1", pipeline randomization removed) and 
deconvolved (bottom right; arestore with 100 iterations) images of the Cas A CCO. See text for details. 



-14 -2 

ergs cm s 



' arcmin ^ for the 
3.45) precludes any de- 



face brightness of ~ 6 x 10 
PL model with F = 1 .5 and nH.22 
tailed analysis of this diffuse emission. 

Source X falls near the edge the EPIC FOV in the 2001 ob- 
servation (Fig. 5b). Although some emission at the source X 
position is discernible, it is not possible to determine whether 
the source is extended or pointlike. Its brightest part in the 
EPIC image is shifted by about 30" southwest from the Chan- 
dra position (the difference could result from a statistical noise 
and background fluctuations). 

Source X was outside the EPIC FOV in the 2005 XMM- 
Newton observation. Because of the low surface brightness 
and large XIS PSF, the source is not seen in the Suzaku data. 
Deeper ACIS or EPIC observations, with the center of the TeV 
source imaged close to the optical axis, are required to under- 
stand the nature of the detected extended X-ray emission. 

In addition to analyzing the X-ray images, we examined 
the available IR and radio data on the region of interest. The 
843 MHz image from the Sydney University Molonglo Sky 
Survey (SUMSS) shows some diffuse emission surrounding 
the very central part of HESSJ1616 (within the innermost 
TeV contour shown in Fig. 5a,f). The radio emission looks 
like a patchy, elongated shell, possibly an unknown SNR (al- 
ternatively, it could be several faint point sources accidently 
aligned in a shell-like structure). The SpitzerlRAC images of 
the same region (Fig. 5e,f) also reveal diffuse emission partly 
coincident with the radio emission. In addition, there are sev- 
eral compact sources seen in the radio and, especially, IR im- 
ages within the central part of the TeV source. The brightest 
source, marked A in Fig. 5, is located near the TeV source 
center At 8 /im, source A is resolved into a complex shell 
with an additional lobe to the west and a number of compact 
point sources suggesting a star-forming region. Source X ap- 
pears to have IR counterpart(s) resolved into several compact 
diffuse sources at 8 /im (Fig. 5e), possibly also star-forming 



regions. Deeper radio observations with higher spatial resolu- 
tion would be most helpful to understand the nature of these 
sources. 

2.3. Pulsar timing 

The data from the CC-mode observation of the RCW 103 
CCO also allow one to analyze the pulsations of the J 1617 
pulsar since the pulsar falls within the FOV. The pulsar is 
clearly seen in the one-dimensional image shown in Figure 6; 
however, the large off-axis angle (w 6') significantly broadens 
the PSF, making it impossible to separate the inner PWN and 
the pulsar 

For the timing analysis, we extracted 4057 photons in the 
2-8 keV band from the 10"-wide segment centered on J1617 
(shown in Fig. 6, top; « 90% of these counts are expected 
to come from the unresolved source, i.e., the pulsar and the 
compact inner PWN). The photon arrival times have been 
transformed to the solar system bary center using the axBary 
tool of ClAO 3.4. We have searched for the pulsed sig- 
nal near the expected radio pulsation frequency and found 
Z?„,^ = 301 (see Buccheri et al. 1983 for details on statis- 
tics) at 1/ = 14.414488 Hz ± l^Hz, consistent with the fre- 
quency expected from the radio ephemeris at the epoch of 
the Chandra observation. The light curve folded at this fre- 
quency reveals a single pulse with a flat minimum (see Fig. 6, 
bottom). The observed 2-8 keV pulsed fraction is 40% ± 4% 
Correcting it for the background contribution (which includes 
the outer PWN), we obtain 44% ±4%. The pulse profiles 
extracted in narrower energy bands (2^ and 4-8 keV) are 
similar to the 2-8 keV pulse profile, both in shape and pulsed 
fraction. 

2.4. Spectral analysis 
2.4.1. Pulsar 



6 



ACIS 



• Kes32 



3 ISUMSS 843 MH 



N. 



MOS1+2 




RCW103 




MOS1+2 


c 















* 



Fig. 5. — X-ray, radio and IR images of the J1617 and HESS J1616 field, a: Combined 0.5-8 keV image from two ACIS observations (2005 June and 2001 
October) with TeV contours (from Aharonian et al. 2006) overplotted. h and c: MOSl-l-2 images in the 0.5-8 keV band from the XMM-Newton observations of 
2001 September and 2005 August, respectively. The cross shows the center of HESS J1616. d: 843 MHz SUMSS image. The rectangles show the region covered 
by the ACIS S2 and S3 chips in the Chandra observation of 2001 October, while the circles show the MOSl+2 FOV in the two XMM-Newton observations. The 
position of J1617 is shown by the diamond, e and/: Spitzer images of the field. 



For our 1/4 subarray observation (ObsID 6684), we ex- 
tracted the pulsar's spectrum (using the ClAO's psextmct 
task) from a small circular aperture of r = 0.5" to minimize 
the contamination from the compact PWN (see Fig. 3, bot- 
tom, and also Fig. Aa,d). No background subtraction was 
attempted because the PWN contribution is expected to be 
negligible in this small aperture (Fig. 3, bottom). The 2853 
counts extracted from the 0.5" radius aperture (52% encir- 
cled energy fraction) in the 1-8 keV range (there are virtually 
no counts below 1 keV) correspond to the aperture-corrected 
source count rate of 96 ±2 counts ks"'. The corresponding 



absorbed flux is /v,. 



2.15 X 10" 



ipsr _ ^.i^ ergs cm s" . With the 

(exposure) frame time of 0.8 s, the expected pile-up fraction 
is very small (« 3% according to PIMMS^). 

We binned the spectrum with minimum of 50 counts per 
bin and fit the absorbed PL model in the 1-8 keV range with 
all the fitting parameters allowed to vary. The quality of the 
fit is excellent (xj, = 1-04 for 49 degrees of freedom [dof]; 
see Fig. 7). The hydrogen column density, «h,22 = 3.45 ±0.14 
(see Fig. 8; the errors here and below are at the 68% confi- 
dence level for a single interesting parameter) is a factor of 

' http://cxc.harvard.edu/toolkit/pimms.jsp 
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Fig. 6. — Results from the CC-mode observation. Top: One-dimensional 
profile showing J1617 and the outer PWN above the horizontal dashed line 
that corresponds to the background level. The two vertical lines show the 
aperture from which the spectrum and the light curve were extracted. Bottom: 
The light curve with 30 phase bins. The dashed line shows the background 
level (including the contribution of the outer PWN). 
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Fig. 7. — Spectrum of the pulsar in observation 6684 fitted with the PL 
model (see Table 3 for the fitting parameters). 
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Fig. 8.— Top: Confidence contours (68%, 90%, and 99%) in the m^-T 
plane for the PL fit to the pulsar spectrum from two ACIS observations (see 
Table 1). The contours are obtained with the PL normalization fitted at each 
point of the grid. Bottom: Confidence contours (68% and 90%) for the PL 
fit to the pulsar spectrum for fixed nH,22 = 3.45. The PL normalization is in 
units of 10"' photons cm~^ s~' keV"' at 1 keV. The dashed curves are the 
loci of constant unabsorbed flux in the 0.5-8 keV band; the flux values near 
the curves are in units of 10"'^ ergs cm~^ s"' . 

1 .5 larger than the average Galactic HI column density in the 
direction of the pulsar (/ = 332.50°, h = -0.28°), nm.ii = 2.26 
(Dickey & Lockman 1990). This is not stirprising since 
the hh deduced from an X-ray spectrum under the assump- 
tion of standard element abundances generally exceeds the 
hhi measured from 21 cm observations by a factor of 1.5- 
3 (e.g., Baumgartner & Mushotzky 2005). The same PL 
fit gives the photon index Fpsr = 1.14 ± 0.06 and the unab- 
sorbed flux F"™*"^ = (3.6 ± 0.1) X 10-'2 ergs cm'^ s"' in the 
0.5-8.0 keV Sand. It corresponds to the isotropic luminos- 
ity Lpsr « 1.8 X 10^X5 ergs s"! (= 1.1 x lO-^^jE), where 
d(,.i = (i/6.5kpc is the distance to the pulsar estimated from 
the dispersion measure DM = 467 cm"^ pc and the Galac- 
tic electron density distribution model by Taylor & Cordes 
(1993)'°. 

We have also extracted the pulsar spectrum from the Ob- 
sID 2957 data taken in the CC mode, using the same region 
(10"-wide segment) as for the light curve (see §2.3), which 
contains 4522 counts in the 1-8 keV band. Note that in this 

The more recent model by Cordes & Lazio (2002) gives a slightly 
larger distance of 6.8 kpc. We caution, however, that distances based on 
such models are not very certain. For instance, HI absorption measurements 
give ~ 20%-30% smaller distances than the Taylor & Cordes model for two 
pulsars within 20° of J1617 (Kaspi et al. 1998). 
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TABLE 3 

PL FITS TO THE PULSAR AND PWN SPECTRA. 



Region «H,22' A/'-s'' T x^/dof ^x.33'-' /-le'' 

Pulsar 3.45 ±0.14 35.4;^'^ 1.14±0.06 51.0/49 17.92±0.07 

Inner PWN [3.45] 4.2;| ^ 1.47±0.24 17.5/19 1.49!S;?^ 940±160 

OuterPWN [3.45] 5.8;|'^ 1.65ifc0.20 12.9/15 1.71ifco'25 1.3ifc0.2 

Note. — The uncertainties are given at the 68% confidence level for one 
interesting parameter. Systematic uncertainties are not included. The nor- 
malization and luminosity of the pulsar are corrected for the finite extraction 
aperture. 

" nil 22 = "h/10'^ cm"2 was fixed in the PWN fits. 
''Spectral flux in units of 10"^ photons cm"- s"' keV"' at 1 keV. 
Isotropic luminosity in units of 10^' ergs s"' in the 0.5-8 keV band, for 
d = 6.5 kpc. 

''Mean unabsorbed intensity in units of 10"'* ergs s"' cm"^ arcsec"^ in the 
0.5-8 keV band. 

case the extraction region is much larger than that used for 
ObsID 6684 because the broad off-axis PSF. The background 
was taken from a nearby (1.4'-2.4' from the pulsar) region lo- 
cated on the same chip (S2). The background contribution is 
« 8% in the 1-8 keV band. The absorbed PL fit gives T ~ 0.85 
and «H,22 — 3.8 (see Fig. 8 for the statistical uncertainties). 
Although the difference between the T values from ObsID 
6684 and ObsID 2957 appears to be statistically significant, 
the CC-mode fit is likely less reliable because of larger un- 
accounted systematic uncertainties (the CTI, gain, and ACIS 
filter contamination corrections are less accurate in the CC 
mode); therefore, below we will use the more reliable spec- 
tral parameters inferred from the ObsID 6684 data. 

2.4.2. PWN 

We used the ObsID 6684 data to extract the PWN spectra 
from two regions. The brightest inner component is extracted 
from the 0.75" < r < 1.25" annulus with the area of 3.14 
arcsec^ (see Fig. 4). In this case we used the pulsar spectrum 
as the background spectrum, after scaling the normalization 
according to the simulated PSF shown in Figure 3 (bottom). 
The fainter outer component is extracted from a polygon re- 
gion (see Fig. 2) with an area of 2600 arcsec^ (excluding the 
r = 9" circle around the pulsar). The background was taken 
from the 40" x 65" box shown in Figure 2 (top). 

The annulus region of the inner PWN component has 1 162 
counts (see Table 2 for details), with w 63% of them coming 
from the pulsar (treated as a background). The corresponding 
absorbed flux is fi„ner = (1.6 ±0.3) x lO"'"* ergs cm"^ s"' in the 
0.5-8 keV band. The outer PWN region contains 959 counts 
(54% coming from the source); its absorbed flux is Fouter = 
(1.7±0.1) X lO"'-' ergs cm"^ s~' in the same band. 

We binned the spectra extracted from the inner and outer re- 
gions with a minimum of 50 counts per spectral bin. We then 
fit the spectra with the absorbed PL model with the hydro- 
gen column density fixed at the best-fit value obtained from 
the PL fit to the pulsar spectrum, nH,22 = 3.45. For the outer 
PWN component, we restrict the energy range to 2-8 keV 
because of the large contribution of the non-uniform back- 
ground" at lower energies (the background contribution de- 
creases from 46% in 0.5-8 keV to 36% in the 2-8 keV band). 
We find that the PL model provides a good description of the 
extracted spectra for both the inner and outer PWN. The pho- 

' ' Investigation of the background spectrum shows that it has a soft com- 
ponent likely originating from a dust-scattering halo around the X-ray bright 
RCW 103. 
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Fig. 9. — Confidence contours (68% and 90%) for the PL fit to the inner 
and outer PWN regions (red and blue contours, respectively; see text for the 
region definitions) for fixed fin, 22 = 3.45. The PL normalization is in units 
of 10"' photons cm"- s"' keV"' at 1 keV. The dashed curves are the loci 
of constant unabsorbed flux in the 0.5-8 keV band; the flux values near the 
curves are in units of 10"'^ ergs cm"- s"' . 

ton indices are Tinner = 1 .47 ± 0.24 and Fouter = 1 -65 ± 20 (see 
Fig. 9). Thus, the PWN spectra are noticeably softer than the 
pulsar spectrum. The combined unabsorbed luminosity from 
the two PWN regions is Lp„n ~ 3.2 x \^^d\^ ergs s~' in the 
0.5-8 keV band (see Table 3 for details). 

3. DISCUSSION 

In our high-resolution Chandra ACIS observation, we dis- 
covered the X-ray PWN around the young, energetic pulsar 
J1617. Below we discuss the pulsar and PWN properties com- 
pare them with those of other young pulsars, and discuss the 
connection between J1617 and the nearby HESS J1616. 

3.1. Pulsar 

The X-ray luminosity of J1617 is ^ 15% lower than that 
reported by BA02 from the XMM-Newton observations. The 
discrepancy is not surprising because the broader PSF of 
XMM-Newton includes both the pulsar and PWN contribu- 
tions. Once these components are added together, the Chan- 
dra and XMM-Newton flux (and spectral slope) measure- 
ments are in agreement within their uncertainties. The 44% 
pulsed fraction measured from the CC-mode data is somewhat 
lower than that reported by BA02. The unsubtracted contribu- 
tion from the inner PWN, resolved in the deconvolved images 
(see Figs. 3 and 4), is « 10%-15%, resulting in the intrinsic 
pulsar's pulsed fraction of w 50%. 

The X-ray efficiency of the J1617 pulsar, 77psr = L^sr/E ~ 
1.1 X IQr^dl^ in the 0.5-8 keV band, is typical among young 
pulsars (see Fig. 10, top). However, the spectrum of J1617 
is surprisingly hard, Fpsi- « 1.1. Very few pulsars (J1509-58, 
J 1420-6048, J 1846-0258) possibly exhibit similai'ly hard X- 
ray spectra (KP08 and references therein), but the accuracies 
of the spectral slope measurements in these pulsars are lower 
than in J1617. The fact that the energy spectrum of the pulsar 
magneto spheric X-ray emission can be so hard has important 
implications for the pulsar emission models. For instance, 
the models that interpret the non-thermal X-ray emission as 
synchrotron emission from the particles produced in a pair 
cascade predict 1.5 < F < 2 (e.g., Cheng et al. 1998; Crusius- 
Watzel et al. 2001), i.e. softer than the spectrum we measured. 
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On the other hand, the so-called "full polar cap cascade" mod- 
els, in which the X-rays are produced by the resonant in- 
verse Compton scattering of the thermal emission from the 
hot NS surface (Zhang & Harding 2000), can explain harder 
spectra, with Tpsr ~ 1 (Fang & Zhang 2006). Hard spectra 
(F ~ 2/3) are possible in the models with curvature-radiation- 
induced cascades in which curvature radiation may dominate 
synchrotron radiation even below ~ 10 keV (Rudak & Dyks 
1999), although this emission mechanism is more commonly 
considered for higher energies. Finally, a hard spectrum could 
be produced via the resonant inverse Compton scattering of 
the radio photons if this process contributes substantially to 
the X-ray band (Petrova 2008). Measurements at higher ener- 
gies with Integral and RXTE suggest a mild softening of the 
pulsar spectrum, by AF « 0.2-0.3 (see Fig. 5 in Landi et al. 
2007). Such behavior is opposite to the spectral hardening ob- 
served in magnetars, which is attributed to comptonization of 
the thermal emission from the NS surface by the relativistic 
electrons in the NS magnetosphere (Fernandez & Thompson 
2007; Baring & Harding 2008). 

The lack of thermal component in the detected pulsar emis- 
sion might be explained by its intrinsic faintness relative 
to the strong magnetospheric emission in this young pulsar. 
However, it is hard to constrain its contribution because, at 
«H 22 ~ 3.5, the soft thermal X-rays are strongly absorbed by 
the ISM. 

3.2. PWN 

Unlike the pulsar, the X-ray efficiency of the J1617 PWN 

(compact + extended), rypwn = Lpv/n/E = 2 x 10~^t/g5 in the 
0.5-8 keV band, is noticeably lower then those of other 
young PWNe (Fig. 10, middle), except for J1357-6429 (Za- 
vUn 2007). Also, the ratio of the PWN to pulsar luminosities, 
Lpvm/LpsT ~ 0.18 (independent of the poorly known distance), 
is among the lowest (Fig. 10, bottom), and it is significantly 
lower than the average Lpwn/ipsr ^ 4 reported by Kargaltsev 
et al. (2007a) and KP08. The actual PWN efficiency and the 
Lpv/n/LpsT value could be higher for J1617 (and J1357-6429) 
if the compact PWN component is not completely resolved 
from the pulsar. An upper limit on the luminosity of the unre- 
solved PWN core follows from the timing analysis (see §2.3 
and §3.1). The luminosity of the unresolved core cannot ex- 
ceed « 50% of the pulsar luminosity (assuming that the pulsar 
emission is 100% pulsed), which still leaves the J1617 PWN 
underluminous compared to majority of PWNe (see Fig. 10 
and the discussion below). 

KP08 found that the X-ray efficiencies of PWNe show a 
very large scatter, rj^^n ^ 10"^-10~' . The cause of such a large 
scatter is still not understood. One could attribute the scat- 
ter to environmental differences (e.g., ISM pressure and den- 
sity); however, the surprisingly good correlation between the 
PWN luminosities and non-thermal pulsar luminosities (with 
a much smaller scatter, 0.1 < Lp^n/^psr ^ 10; KP08) calls for 
a different explanation because Lpsr does not depend on the 
environment. 

The deconvolved Chandra images of J1617 suggest a com- 
pact structure at a distance of ^ 1" (= 9.7 x lO'^'fifi.s cm) from 
the pulsar (see Fig. 4), which could be a torus beyond the ter- 
mination shock in the pulsar wind. Assuming that the ram 
pressure caused by the pulsar motion is lower than the ambi- 
ent pressure, we can estimate the latter as 

Pamb ~ Efai47rcr^J-' (1) 
where is the termination shock radius in the equatorial 
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Fig. 10. — Top: Pulsar luminosity vs. spin-down power £. Middle: PWN 
luminosity vs. E. Bottom: Pulsar luminosity vs. PWN luminosity. The top 
and middle panels show the lines of constant X-ray efficiency, while the 
dashed line in the bottom panel corresponds to Lpsr = Lpwn- (Based on the 
data pubUshed by KP08.) 

plane, and the factor /j^ takes into account possible anisotropy 
of the pulsar wind. For J1617, this estimate gives p-^^b ~ 
1.7 X 10"Vf2(rv/5 X 10'^' cm)"^ ergs cm"^ It is unhkely that 
the X-ray PWN is more compact than it follows from our im- 
age analysis (which would imply that the true PWN lumi- 
nosity and its ratio to Lpsr are higher than estimated above) 
because the inferred pressure is already quite high (cf. Table 
2 in Kargaltsev et al. 2007a). Therefore, we suggest that in 
some PWNe, including J 1617, the radiative efficiency of the 
shocked pulsar wind is well below the average. The low PWN 
luminosity cannot be attributed to inefficient pair production 
in the pulsar magnetosphere because the non-thermal pulsar 
luminosity is close to its typical value. However, the faintness 
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of the PWN could be attributed to inefficient particle accelera- 
tion on the way to (or at) the termination shock. The physical 
mechanism responsible for particle acceleration is currently 
unknown, and the processes occurring upstream of the shock 
are poorly understood (Arons 2007; Kirk et al. 2007). We 
can only speculate that these processes may depend of the 
angle between the pulsar magnetic and rotation axes (e.g., a 
larger misaUgnment of the axes may lead to faster conversion 
of the magnetic field energy to the particle energy via more 
efficient recormection). They may also depend on the termi- 
nation shock radius, which, in turn, depends on the ambient 
pressure. For instance, if is small, there may not be enough 
time for the dissipative processes to convert the magnetic field 
energy to the particle energy (e.g., Arons 2008). Yet, even 
in this case, most of the pulsar's rotational energy is being 
lost as a wind, but the wind may have a higher magnetization 
parameter and lower particle energies (see, e.g., Arons 2007 
for a pulsar wind theory review) than in the case of X-ray 
bright PWNe. A higher-than-typical ambient pressure (hence 
a small termination shock radius) is plausible because the 
young pulsar should be within the host SNR interior (which 
remains undetected because of the strong X-ray absorption). 
The wind still should emit synchrotron radiation downstream 
of the termination shock but possibly at lower frequencies as 
the characteristic synchrotron frequency is proportional to the 
electron Lorentz factor squared. Unfortunately, observing at 
lower X-ray energies or in the optical is impossible because 
of the strong ISM absorption, but sensitive far-IR and radio 
observations would be useful. 

3.3. Connection to HESS J1616-508 

J 16 17 is located - 10' east of the TeV source HESS Jl 616- 
508 (see Fig. 5). Despite the large offset, the yoimg, ener- 
getic pulsar is still within the TeV source extent, and it has 
been considered as a plausible candidate for powering the TeV 
emission (Aharonian et al. 2006; Landi et al. 2007), mainly 
based on energetics arguments (Lxev ~ lO"^^). Extended TeV 
sources have recently been discovered in the vicinity of a few 
(w 15) young pulsars, leading to the suggestion that the TeV 
emission comes from relic PWNe crushed by the asymmetric 
SNR reverse shock that arrived to the pulsar location a few 
kyrs ago (e.g., de Jager 2006; Kargaltsev et al. 2007b; Kar- 
galtsev & Pavlov 2007). The original model developed by 
Blondin et al. (2001) suggests that the synchrotron emitting 
electrons from the pulsar wind can be swept up by the asym- 
metric reverse shock to one side of the pulsar. The resulting 
"offset PWNe" will contain aged electrons that still can pro- 
duce observable radio synchrotron emission as well as high- 
energy 7-ray emission via the inverse Compton scattering of 
background photons by relativistic electrons. 

The high-resolution ACIS observations of J1617 do not re- 
veal a preferential extension of the J1617 PWN towards the 
TeV source (rather there is a hint of extension in the opposite 
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direction; see §2.1 and Fig. 1), such as observed in some other 
TeV sources neighboring young pulsars (Kargaltsev & Pavlov 
2007; Pavlov et al. 2008). Therefore, we conclude that there 
is no convincing evidence linking the TeV source to J1617. 

The HESSJ1616 field has been observed with several X- 
ray missions. The XMM-Newton and Suzaku images do 
not reveal any X-ray sources within the TeV bright part of 
HESS J1616, making it one of the best examples of a "dark 
accelerator" (Matsumoto et al. 2007). Yet, the examination 
of the 30 ks ACIS image covering the central part of the 
HESS J1616 field, reveals the possibly extended source X (see 
Fig. 5 and §2.2 ). This source, however, has a very low X-ray 
flux resulting in a surprisingly small Lx /^xev = (1 ~2) x 10~^ 
(cf. Table 2 in Kargaltsev et al. 2007b). Furthermore, the radio 
and IR images shown in Figure 5 reveal some diffuse emission 
and pointlike sources within the central part of HESS J1616. 
The diffuse emission resembles a shell, while one of the bright 
IR pointlike sources coincides with the center of the TeV 
source. Therefore, we speculate that HESS Jl 61 6 may con- 
sist of two (or more) unresolved sources, one of which stiU 
might be associated with J1716, while the other(s) could be 
associated with the extended radio/IR emission (possibly an 
SNR or a star-forming region). Deeper radio and X-ray ob- 
servations or TeV observations with higher resolution can test 
this hypothesis. 

4. SUMMARY 

We have discovered a surprisingly underluminous PWN 
around the young, energetic pulsar JI617. The PWN consists 
of a compact and extended components of comparable lumi- 
nosities. The faintness of the PWN could be attributed to the 
intrinsically low radiative efficiency of the pulsar wind. We 
have also obtained an accurate measurement of the pulsar's 
X-ray spectrum, which fits an absorbed PL with Fpsr 1.1, 
harder than predicted by the models of synchrotron emis- 
sion from relativistic electrons/positrons produced in the pair- 
cascade. 

J1617 has been considered as a plausible counterpart to 
file offset TeV source HESS J1616; however, tiie X-ray PWN 
does not show any asynnmetry towards HESS J1616. On the 
other hand, the archival X-ray, IR, and radio data that cover 
the central region of HESS J 1616 reveal some diffuse emis- 
sion, allowing for altemative interpretations of the TeV emis- 
sion. 
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